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Abstract. The spectrum of the Laplacian of successive Barycentric 
subdivisions of a graph converges exponentially fast to a limit which 
only depends on the clique number of the initial graph and not on the 
graph itself. Announced in [40]), the proof uses now an explicit linear 
operator mapping the clique vector of a graph to the clique vector of the 
Barycentric refinement. The eigenvectors of its transpose produce inte¬ 
gral geometric invariants for which Euler characteristic is one example. 


1. Notations 

Given a finite simple graph G = {V, E) with vertex set V and edge set E, the 
Barycentric refinement Gi = (Vi,Ei) is the graph for which Vi consists 
of all nonempty complete subgraphs of G and where Ei consists of all un¬ 
ordered distinct pairs in Vi, for which one is a subgraph of the other. Denote 
by Gjn the successive Barycentric rehnements of G assuming Go = G. If 
^0 ^ ^ -^n are the eigenvalues of the Kirchhoff Laplacian L = B — Aoi 

G, where B is the diagonal degree matrix and A the adjacency matrix 
of G, define the spectral function Eq{x) = where [t] is the largest in¬ 
teger smaller or equal to t. So, Eg{0) = 0 and Eg{1) is the largest eigenvalue 
of G. The norm of F satisfies IILgHi = ||A||i/|V| = tr(L)/|V| = 2|E|/|V| 
by the Euler handshaking lemma telling that 2|E| is the sum of the ver¬ 
tex degrees. In other words, ||E||i = d{\V\ — 1), where d is the graph 
density d = \E\/B{\V\,2) with Binomial coefficient B{-,-)- The density d 
gives the fraction of occupied graphs in the completed graph with vertex set 
V. If G, H are two sub graphs of some graph with n vertices, define the 
graph distance d{G, H) as the minimal number of edges which need to 
be modified to get from G to H. If L, K are the Laplacians of G, H, then 
j l-bij “ Kij\ ^ 4d(G, R) because each edge (i,j) affects the four matrix 
entries La, Ljj of the Laplacian only: adding or removing that edge 

changes each of the 4 entries by 1. The Lidskii-Last inequality assures 
IIai — A||i < Y17j=i 1^ “ [SSI US] for any two symmetric n x n matrices 
A, B with eigenvalues ai < 0:2 < • • • < and Pi < j32 < ■ ■ ■ < Pn- For two 
subgraphs G,H of a common graph with n vertices, and Laplacians L, H, 
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the inequality gives ||A —//||i < 4:d{G, H) so that \\Fg' — Fh'\\ < 4:d{G, H)/n 
if G', H' are the graphs with edge set of G and vertex set of the host graph 
having n(G, if) vertices. Therefore \\Fg — Fh\\i < 4:d{G, H)/n{G, Ff) holds, 
where n(G, H) is the minimum of |F(G)) and \ V{F[)\ assuming both are in a 
common host graph. The distance of spectral functions can so estimated 
by the graph distance. Obviously, if k disjoint copies of the same graph H 
form a larger graph G, then F{G) = F{H). Given a cover Uj of G, let H be 
graph generated by the set of vertices which are only in one of the set Uj . Let 
K be the graph generated by the complement of H. Now, d{G,H) < 4|iL| 
and d{Uj, Hj) < 4|iL|. f is the number of complete subgraphs iLfc+i of 
G, the clique number of G is defined to he k if Vk = 0 and Vk-i > 0. A tree 
for example has clique number 2 as it does not contain triangles. Denote by 
Qd the class of graphs with clique number d + 1. The class contains Kd+i 
as well as any of its Barycentric refinements. The clique data of G is the 
vector V = (uq, ui,...,), where Vk counts the number of K^+i subgraphs of 
G. We have vq = |D|,ui = \E\ and V 2 counts the number of triangles in 
G. The clique data define the Euler polynomial edx) = and 

the Euler characteristic x(G) = eG'(—1). The polynomial degree of e(x) 
is d if d + 1 is the clique number. There is a linear Barycentric opera¬ 
tor A which maps the clique data of G to the clique data of Gi. It is an 
upper triangular linear operator on 1“^ with diagonal entries Akk = k\. It 
also maps the Euler polynomial of G linearly to the Euler polynomial of Gi. 
The unit sphere S{x) of a vertex x G E in G is the graph generated by 
all vertices connected to x. The dimension of G is defined as dim(0) = —1 
and dim(G) = 1 + Y.v(iV dim(S'(x))/uo, where S{x) is the unit sphere. A 
graph has uniform dimension d if every unit sphere has uniform dimen¬ 
sion d — 1. The empty graph has uniform dimension —1. Given a graph G 
with uniform dimension d, the interior is the graph generated by the set of 
vertices for which every unit sphere S{x) is a (d — l)-sphere. A d-sphere 
is inductively defined to be a graph of uniform dimension d for which every 
S{x) is a (d — l)-sphere and for which removing one vertex renders the graph 
contractible. This Evako sphere definition starts with the assumption that 
the (—l)-sphere is the empty graph. Contractibility for graphs is induc¬ 
tively defined as the property that there exists x G E for which both S{x) 
and the graph generated by E \ {x} are both contractible, starting with the 
assumption that Ki is contractible. The boundary 5G of a graph G is the 
graph generated by the subset of vertices in G for which the unit sphere S{x) 
is not a sphere. Also the next definitions are inductive: a d-graph is a graph 
for which every unit sphere is a (d— l)-sphere; a d-graph with boundary 
is a graph for which every unit sphere is a (d— l)-sphere or (d— l)-ball; a d- 
ball is a d-graph with boundary for which the boundary is a (d— l)-sphere. 
Let Wk{Gm) denote the number of Kk+i subgraphs in the boundary 5Gm- 
For G = Kd+i, the boundary of Gm is a (d — l)-sphere for m > 1 and 5Gm 
contains the vertices for which the unit sphere in Gm has Euler character¬ 
istic 1. Barycentric refinements honor both the class of d-graphs as well as 
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the class of d-graphs with boundary. Starting with G = Kd+i which itself is 
neither a d-graph, nor a d-graph with boundary, the Barycentric refinements 
Gm are all d-balls for every m > 1. While the spectral functions F are well 
suited to describe limits in L^([0,1]), one can also look at the integrated 
density of states a monotone [0, l]-valued function on [0,oo) which 
is also called spectral distribution function or von Neumann trace. 
It defines the density of states {F~^)' which is a probability measure on 
[0, oo) also called Plancherel measure, analogue to the cumulative distri¬ 
bution functions defining the law of the random variable which in absolutely 
continuous case is the probability density function. Point-wise convergence 
of F implies point-wise convergence of F~^ and so weak-* convergence of 
the density of states. 

2. The theorem 

Theorem 1 (Central limit theorem for Barycentric subdivision). The func¬ 
tions Fg^{x) converge in T^([0,1]) to a function F{x) which only depends 
on the clique number of G. The density of states converges to a measure p, 
which only depends on the clique number. 

We first prove a lemma which is interesting by itself. It allows to compute 
explicitly the clique vector of the subdivision Gi from the clique vector of 
G. 

Lemma 2. There is an upper triangular matrix A such that v{Gi) = Av{G) 
for all finite simple graphs G. 

Proof. When subdividing a subgraph ii'fc+i h splits into A^^ = {k + 1)! 
smaller graphs. This is the diagonal element of A. Additionally, for 

m > k, every K^+i subgraph produces A^m subgraphs isomorphic to A'fc+i) 
which is the number of interior subgraphs of the Barycentric subdi¬ 

vision of Kd+i. These ATm+i subgraphs in the interior correspond to Km +2 
subgraphs on the boundary. This means that Af^^ is the number of Km +2 
subgraphs of the boundary of the Barycentric refinement of Since 

this boundary is of smaller dimension, we can use the already computed 
part of A to determine A^m- To construct A, we build up the columns re¬ 
cursively, starting to the left with ei. Let B{n,k) = nl/{kl{n — k)\). If n 
columns of A have been constructed, we apply the upper n x n part of A 
to the vector v = [B{n -|- 1,1),..., B{n -|- 1, n)]'^ in order to get the clique 
data of the (n — l)-sphere 5{Kn+i)i. These numbers encode the number of 
interior complete subgraphs of the n-ball (iLn+i)i- If the resulting vector is 
[rci,..., Wn]'^, take [1, rci,..., Wn-i, {n -|- 1)!, 0,0,... as the new column. 
The Mathematica code below implements this procedure. □ 

The Barycentric operator A has an inverse A~^ which is a bounded 
compact operator on f^(N). Below we give the bootstrap procedure to 
compute A, adding more and more columns, using the already computed 
A to determine the next column. The eigenvalues A of A are included in 
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the spectrum a{A) = {!!, 2!, 3!,... }. Any eigenvector / of can lead to 
invariants X{G) = {f,v{G)) which correspond to valuations in the con¬ 
tinuum as XX(G) = {f,v) = {A^f,v) = {f,Av) = X{Gi). Of particular 
interest is the Euler characteristic eigenvector / = [1, —1,1, —1, ....]^ to 
the eigenvalue A = 1 verifying that x(G') is an invariant under Barycentric 
subdivision. An other invariant is the eigenvectors / = [0,..., 0, —2, d + 1]^ 
to the eigenvalue d\ which measures the “boundary volume” of a d-graph 
with boundary. For d-graphs, it remains zero under refinement. 

Here is the proof of the theorem: 

Proof, (i) For G G Gd^ there are constants > 0 such that Vk{Go){d+ 
If-T/Cd < Vk{Gm) > Ufc(Go)Crf((d+l)!r andt(;fc(Go)(d!r/Crf < Wk{Gm) > 
WkiGo)Gd{d\)'^ for every 0 < k < d. Proof. This follows from the lemma 
and Perron-Probenius as the matrix A is explicit. When restricting to Qd, 
it is a (d -|- 1) X (d -|- l)-matrix with maximal eigenvalue (d -|- 1)!, when re¬ 
stricted to Gd-i its maximal eigenvalue is d!. Diagonalizing Ad = Sf^BdSd 
we could find explicit bounds Gd- 

(ii) For every G = K^, the sequence converges. Proof: There 
are (d-|- 1)! subgraphs of Gm=i isomorphic to Gm-, forming a cover of Gm+i- 
They intersect in a lower dimensional graph. Since the number of ver¬ 
tices of this intersection grows with an upper bound Gddl^ in each Gm we 
have d(Gm+i,UjG'm^) < {d + l)!C'd(d!)”*. This shows \\Fg^ - Fb^+il|i < 
Gd{d + l)!(d!)”^/((d-|- 1)!)™ = Gd/{d + 1)™“^. We have a Cauchy sequence 
in T^([0,1]) and so a limit in that Banach space. 

(ii) Barycentric refinement is a contraction on Gd in nn adapted 
metric. Proof. Every G & Gd can be written as a union of several Hm ~ 
Kd+i subgraphs and a rest graph L in Gk with k < d. Then Fq^ and 
Fh^ have the same limit because Fh^ = Fb^-i- The Barycentric evolution 
of the boundary of refinements of Kd+i as well as L grows exponentially 
smaller. Given two graphs A, B G Gd, then d{Am, Bm) < Cd(d!)”^ so that 
I|Ta^-Fb„||i < 4Crf/(d-|-l)™'. Letmobesolargethatc = C"/(d-|-l)™'° < 1. 
Define a new distance d'{A,B) = ^i^k, Bk), so that d'{Ai,Bi) < 

cd'{A,B). Apply the Banach fixed point theorem. 

(vi) There is uniform convergence of Fbm ®n compact intervals 

of (0,1). Since each Fq^ is a monotone function in L^([0,1]), the ex¬ 
ponential convergence on compact subsets of (0,1) follows from exponen¬ 
tial convergence. (This is a general real analysis fact [IS]). In di¬ 
mension d, one has ||Tb,„||i —)• (d -|- 1)! exponentially fast. Indeed, as 
the number of boundary simplices grows like (d!)™ and the number of 
interior simplices grows like ((d -|- I)!)™', the convergence is of the order 
l/(d-|-l)'". By the Courant-Fischer mini-max principle Fb(l) = Xn = 
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viiax.{v,Lv)/{v,v) > max(La;a;) = max(deg(x)) grows indefinitely, so that 
the L°° convergence can not be extended to 1 ]. □ 

It follows that the density of states of G converges “in law” to a universal 
density of states which only depends on the clique number class of G, hence 
the name “central limit theorem”. The analogy is to think of G or its Lapla- 
cian L as the random variable and the spectrum <t{L) of the Laplacian L as 
the analog of the probability density and of the Barycentric refinement oper¬ 
ation as the analogue of adding and normalizing two independent identically 
distributed random variables. For d = 1, where the graph is triangle-free 
and contains at least one edge, we know everything: 

Proposition 3. For d = 1, the limiting function is Ti(x) = 4 sin^( 7 rx/ 2 ). 

Proof. As the limiting distribution is universal, we can compute it for G = 
Cm, where Gm = G 4 . 2 ™. As the spectrum of Cn is the set {Asiv?{'Kk/n) \ k = 
1 ,..., n}, the limit is Fi. □ 

For d = 1, the limiting spectral function is related to the Julia set of 
the quadratic map 2 ; —)• z(4 — z) which is conjugated to z —?■ z^ — 2 
(c = —2 at bottom tail of Mandelbrot set) or 2 ; —)• 42;(1 — z) which is the 
Ulam interval map conjugated to tent map, or z —?■ 42;^ — 1 which is a 

Chebyshev polynomial. 


3. Remarks 

1) For d = 2 already, we expect spectral gaps. A first large one is observed 

at X = 1/2. For G 4 with G = K 3 , we see a jump at 0.5 of 2.002, Starting 
with G = ATs, the graph G 2 has 25 vertices. Its Laplacian has eigenval¬ 
ues for which A 13 — A 12 = 2.0647... ~ 2 being already close to the gap. 
The eigenvalues can be simplified to be roots of polynomials of degree 4 for 
which explicit radical expressions exist allowing to estimate. We also know 
||Fg 3 - F 2 II 1 < 1/3^ = 8/18 < 1/2. While we know that Fg™ - F 2 

converges pointwise and uniformly on each compact interval we don’t have 
uniform constant bounds on each interval. 

2) Instead of the Laplacian L = B—A, we could take the adjacency matrix 
A. We could also take the multiplicative Laplacian L' = AB~^ or the to 
L' isospectral selfadjoint random walk Laplacian L” = B^^"^AB^^"^ . The 
corresponding spectral functions always converge exponentially fast, but to 
different limiting functions. For the adjacency matrix for example, the spec¬ 
tral gaps appear much smaller. 

3) The convergence also works for the Dirac operator D = d + d*, where d 
is the exterior derivative or the Hodge Laplacian L = D^. For the Dirac 
operator D, for which the density of states is supported on (— 00 , 00 ). For the 
Hodge Laplacian on [0, 00 ). The Hodge Laplacian factors into different form 
sectors and the spectral functions of each form Laplacian converge. 
The blocks which make up . The scalar Laplacian Fq = d*d agrees with 
the combinatorial Laplacian = Kirchhoff Laplacian L = B — A discussed 
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above. 

4) Any Barycentric refinement preserves the Euler characteristic x(G) = 

We know that Gn is homotopic and even homeomorphic 
to G. G and H are Ivashchenko homotopic [281 [8] if one get H from 
G by homotopy transformation steps done by adding a vertex, connecting 
it to a contractible subgraph, or removing one with a contractible sphere. 
A topology O on V of G is a graph topology if there is a sub-base B 
of O consisting of contractible subgraphs such that the intersection of any 
A,B & B satisfying dim(AnB) > min(dim(A), dim(B)) is contractible, and 
every edge is contained in some B ^ B. We ask the ^ of B to be homotopic 
to G, where the nerve graph Q = {B,£) has edges £ consisting of all pairs 
{A,B) ^ B X B for which the dimension assumption is satisfied. A map 
between two graphs equipped with graph topologies is continuous, if it is 
a graph homomorphism of the nerve graphs such that dim((/>(A)) < dim(A) 
for every A G B. If has a continuous inverse it is a graph homeomor- 
phism. 

5) We know dim(Gi) > dim(G) with equality for d-graphs, graphs for which 
unit spheres are spheres. For a visualization of the discrepancy on random 
Erdos-Renyi graphs, see [l3]. We also know that dim(Gm) converges mono- 
tonically to the dimension of the largest complete subgraph because the 
highest dimensional simplices dominate exponentially. 

6 ) If G is a d-graph, then for m > 1, all Gm as well as any finite intersections 
of unit spheres in Gm are all Eulerian so that we can define a geodesic 
flow on Gm or a billiard on Barycentric refinements of the ball {Kk)m- For 
non-Eulerian graphs like an icosahedron, a light propagation on the vertex 
set is not defined without breaking some symmetry. Also, for any graph and 
m > 0, the chromatic number of Gm is the clique number of G. Indeed, 
for m > 1 the dimension of the complete subgraph graph x of Gm can be 
taken as the color of the vertex x in Gm+i- Since the dimension takes values 
in {0,..., d} the chromatic number of G agrees with the clique number d-|-1 
of G. 

7) If do < • • • < dn is the degree sequence of G, define the degree 
function Hc{x) = dj^.^] in the same way as the spectral function. Since 
the degrees are the diagonal elements of L, the integrated degree func¬ 
tion H{x) = Jq Hcit) dt and the integrated spectral function F{x) = 
fg FG{t) dt agree at x = 0 and x = 1. By the Schnr inequality, we have 
H{x) > F{x). The degree function also converges. 

8) Barycentric refinements are usually defined for realizations of a simpli- 
cial complex in Euclidean space. The Barycentric subdivision of a graph is 
not the same than what is called a simplex graph [3] as the later contains 
the empty graph as a node. It also is not the same than the clique graph 
|24j as the later has as vertices the maximal complete subgraphs, connecting 
them if they have a non empty intersection. The Barycentric refinement is 
the graph product of G with Ki [33]. The product G x H has as vertex 
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set the set of all pairs {x,y), where x is a complete subgraph x of G and a 
complete subgraph y of H. Two such vertices (x, y) and {u, v) are connected 
by an edge, if either x C u and y d v oi u <Z x and y C v. 

9) The limit of Barycentric refinements is “flat” in the following sense: 
the graph curvature K{x) = 1 — Vq{x)/2 + Vi(x)/3 — • • • with V{x) = 
(Vb(3:), ti(x),...) being the clique data of the unit sphere S{x) is defined 
for all finite simple graphs [36] satisfying the Gauss-Bonnet-Chern the¬ 
orem adds up to the Euler characteristic of the original graph. Because the 
total curvature is 1, and as the graph gets larger, the curvature averaged 
over some subgraph Gm-k of Gm goes to zero, while the individual curva¬ 
tures can grow indefinitely. For d = 2 for example, where the curvature 
is K{x) = 1 — d(x)/6, there are vertices with very negative curvature, but 
averaging this over a smaller patch gives zero. In general, if we look at 
subgraphs of G which are wheel graphs as 2-dimensional sections, also the 
sectional curvatures become unbounded at some points but averages of 
sectional curvatures over two dimensional surfaces obtained by Barycentric 
refinements of a wheel graph go to zero. The limiting holographic object 
looks like Euclidean space. There is even “rotational symmetry”, not 
everywhere, but with centers at a dense subset of the continuum. 

10) The graph theoretical subdivision definition uses the point of view that 
complete subgraphs of a graph are treated as points. When this identifi¬ 
cation is taken seriously and iterated, we get a holographic Barycentric 
refinement sequence which is already realized in the graph itself. This 
is familiar when building up number systems: if the set of integers G = Z 
is considered to be a graph with adjacent integers connected, the Barycen¬ 
tric refinements Gn contain all dyadic rationals /c/2"'. Modulo 1, this 
is a Priifer group P 2 which has as the Pontryagin dual the compact 
topological group D 2 of dyadic integers, a subring of the held Q 2 of 2-adic 
numbers. In the case d = 1, there is a limiting random operator on the 
group of dyadic integers. 

11 ) There is an analogy between (M, T, Z, x —>■ x -|- a mod 1 , x —)• 2 x mod 1 ) 

and (Q 2 , IP 2 ) Ih’ 2 ) a: —>■ x -|- 1 , u), where x —>■ x -t- 1 is the addition on the com¬ 
pact topological group. This translation on D 2 is called adding machine 
[20] and a is shift which is a return map on half of D 2 . There is a natural 
way to get the group D 2 through ergodic theory as it is the unique fixed 
point of 2:1 integral extensions in the class of all dynamical system. It 
is called also the von Nenmann-Kakutani system and usually written 
as an interval map on [0,1]. The ergodic theory of systems with discrete 
spectrum is completely understood and the von-Neumann-Kakutani 

system belongs so to a class of systems, where we can solve the dynamical 
log-problem: given any x, y and e > — find n such that d{T'^x, y) < e which 
is essential everywhere in dynamics, from prediction of events up to finding 
solutions to Diophantine equations. Such systems are typically uniquely 
ergodic and so by spectral theory naturally conjugated to a group transla¬ 
tion on a compact topological group, the dual group of the eigenvalues of the 
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Koopman operator on the unit circle in the complex plane. There is an 
important difference between the real and 2-adic story: while in both cases, 
the chaotic scaling systems are isomorphic Bernoulli systems on compact 
Abelian topological groups, the group translation on the 2-adic group is 
naturally unique and quantized, while on the circle, there are many group 
translations x ^ x + a. In the real picture, there is a continuum of nat¬ 
ural translations, in the dyadic picture, there is a smallest translation. 
The picture is already naturally quantum. Egyptian dyadic fractions, music 
notation or the transition from Fourier to wavelet theory can be seen as a 
move towards dyadic models. See [62] for a plethora of other places. 

12) For a d-graph G with d > 1, the renormalization limit of G is still 
unidentified. The limiting operator is likely a random operator on a com¬ 
pact topological group. As the group of dyadic integers, it is likely also a 
solenoid, an inverse limit of an inverse system of topological groups. There 
is a group acting on it producing the operator in a crossed product G*- 
algebra. This is at least the picture in one dimensions. 

13) As Fc^ converges in we have convergence of the integrated den¬ 
sity of states Fq^ in L^[0,cx)) and so weak-* convergence of the deriv¬ 
ative, the density of states, in the Banach space of measures. In the 
case d = 1, the density of states is f{x) = (x(4 — x ))~^/‘^/tt supported on 
[0,4]. The integrated density of states is F~^{x) = (2/7r) arcsin(-y/x/2). The 
measure ^ = /(x)l[o, 4 ]dx is the equilibrium measure on the Julia set of 
T{z) = Az — z^ whose dynamic is conjugated to z —)■ 2 :^ — 2 in the Mandelbrot 
picture or the Ulam map z —?• 4x(l — x) for maps on the interval [0,1]. The 
spectral function F satisfies T{F{x)) = F{2x). The measure ^ maximizes 
metric entropy and equals it to topological entropy log(2) of T. 

14) In the case d = 2, the numbers Vk were already known |60l [5] as 
V 2 {Gm) = S'” and vq is the sum of vo(m — l),vi(m — l),V 2 (m — 1) and 
Vo — vi + V 2 = 2 leading to formulas like vo{Gm) = 1 — 3(2’”“^ + 2"^ + 
2'”-i3'”-i), vi{Gm) = 3(-2'”-i + 2^” + 2'”-i3'”). As the lemma shows, in 
general, the clique data of Gi are a linear image of the clique data of G with 
a linear map A independent of G. Since A has a compact inverse, we can 
look at an eigenbasis of A and could write down explicit formulas for the 
number of vertices of Gn, if the initial clique data of G are known. 

15) If dd is a subgraph of G, then each rehnement Hm is a subgraph of Gm- 
Also, Gk is a subgraph of Gm A k < m and the automorphism group of Gm 
contains the automorphism group of G. The case G = Kn shows that the 
automorphism group can become larger. In general, Aut(Gi) = Aut(G). 
It is only in rare cases like G = Gn that Aut(Gm) can grow indefinitely. 

16) The matrix A mapping the clique data of G to the clique data of Gi is 
a special case of the following: for any graph dd, there is a linear map A^ 
which maps the clique data of any graph G to the clique data of G x FI. 
This operator Ah depends only on dd. While A = Ak^ is invertible on each 
class of finite dimensional graphs. The spectral data v{G) do not determine 
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the graph G as trees already show, the Barycentric refinement operator T 
is invertible on the image G{Q) of the class G of all finite simple graphs. 

17) Barycentric refinements make graphs nicer: the chromatic number is 
the clique number, the graphs are Eulerian. If we identify graphs with their 
Barycentric refinements, we have a holographic picture. There is a nat¬ 
ural geodesic flow on Gi if G is a d-graph as each every unit sphere has a 
natural involution, defined inductively by running the geodesic flow on the 
unit sphere. On a two dimensional sphere for example, the flow is defined 
because every vertex has even degree. This allows to draw straight lines 
on it and to define an antipodal map. For a 3-sphere, because every unit 
sphere has such an antipodal map, we can define straight lines there, in 
turn leading to an antipodal map on such spheres. This can be continued 
by climbing up the dimensions to get a geodesic flow on the graph Gi itself. 

18) The density of states obtained as a limiting spectral density of finite 
dimensional situations is central in the theory of random Schrddinger 
operators L [TTl |53l |50], where the density of states /r can be defined as 
the functional / —)• tr(/(L)) as L is an element in a von Neumann algebra 
with trace, the crossed product of L°°{Q) with an ergodic group action on 
the probability space. The BirkhofF ergodic theorem has been applied 
by Pastur to assure in that theory that the density of states of finite di¬ 
mensional operators defined by orbits starting at some point x converges 
almost surely to the measure fi. The 1-dimensional discrete case d = 1, 
where a single ergodic transformation T on a probability space is given is 
one of the most studied models, in particular the map x ^ x + a mod 1 
leading to almost periodic matrices like the almost Mathieu operator. 
It is the real analog of the limiting operator which is almost periodic over 
the dyadic integers. There is no doubt that also in higher dimensions, the 
limiting model of Barycentric refinement is part of the theory of aperiodic 
almost periodic media IMIIIIIEQIETIEHIEI]. 

4. Open ends 

1 ) Already for d = 2, we do not know what the nature of is. One must 
suspect that there is a renormalization picture, as in one dimensions. It can 
not be a Julia set of a polynomial, as the spectrum becomes unbounded. We 
suspect that for d > 1, the time of the group action is higher dimensional 
and that the limiting Laplacian remains almost periodic. Since we can not 
compute F 2 {Gm) yet for large m due to the fast growth of Gm, we see ex¬ 
perimentally that the function F 2 on L^([0,1]) has a self-similar nature 
in the sense that F 2 { 6 x) ~ (p 2 {F 2 {x)) for some 4>2 and x G [0,1/6]. This 
self-similar nature is present in d = 1 as Fi{2x) = 4’i{Fi{x)), where (pi is 
the quadratic map (pi{x) = x{A — x). Since ^ 2 ( 1 ) is infinite, the map (p 2 
(if it exists at all) must be an unbounded function and in particular can not 
be a polynomial. 
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2) In order to establish an existence of a spectral gap for F 2 , one could 
first establish an exact hnite approximation bound, then estimate the rest. 
For d = 2, the hrst large gap opens at x = 1/2. This could be related to 
the fact that every step doubles the vertex degrees in two dimensions. Here 
are the largest spectral jump values for G 5 if G = K^-. A(0.5) = 2 . 002 .., 
A(0.835366) = 1.67669... A(0.917683) = 2.86249, A(0.972561) = 3.89394, 
A(0.98628) = 6.93379 A(0.995427) = 7.96794 and A(0.998476) = 14.9767. 
To investigate gaps, we would need a Lidskii type estimate for the smaller 
half of the eigenvalues. 

3) We see experimentally that for d = 2, the sequence appears mono- 

tonically increasing Fq^ < at least for the small m in which we can 

compute it. In the case d = 1, the monotonicity is explicit. It fails for 
smaller m and d = 3 but could hold for any d if m is large enough. If we 
had more details about the convergence in the middle of the spectrum, we 
could attack the problem of verifying that spectral gaps exist. 

4) Unlike the Feigenbaum renormalization in one dimensional dynamics 

[laETKIS], where a hyperbolic attractor with stable and unstable manifold 
exists, the Barycentric renormalization is a contraction and the convergence 
proof is more elementary. More so than the central limit theorem in prob¬ 
ability theory, where the renormalization map is A —)• Xi + X 2 , where Xi 
are independent random variables on a probability space H with the same 
distribution than X and X = (X — E[A])/cr(A). That random variable 
renormalization is not a uniform contraction in L^(H,P). For d = 1 there is 
a limiting almost periodic operator where the compact topological group of 
dyadic integers becomes visible in the Schrodinger case. We expect also 
in higher dimensions, a random limiting operator and that the renormaliza¬ 
tion map is robust in the sense that one can for example shift the energy E 
and have a deformed attractor. In the one-dimensional Schroedinger case, 
where one deals with Jacobi matrices, one gets so to spectra on Julia sets 
Je of the quadratic map T{x) = + E [33 ED ED ED ESj- Also in higher 

dimensions, we expect that a Schrodinger renormalization picture reveals 
the underlying topological group if an energy parameter is used to modify 
the renomalization. 

5) In the case d = 1, the roots of the Dirac zeta function (/(s) = X]a>o 
dehned by the positive eigenvalues A of the Dirac operator D = d+d* [381137] 
of Gm, converges to a subset of the line Re(s) = 1. What happens in the 
case d = 21 Already for d = 1, the convergence of the roots is slow of the 
order log(log(uo(G)) as it is initially proportional to m and slowing down 
exponentially when approaching the line. Any experimental investigations 
in d = 2 would be difficult as no explicit formulas for the eigenvalues exist. 
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6 ) It would be nice to know more about the linear operators Ah which be¬ 
long to the map G ^ G x H on Q. Their inverse is compact. In the case 
H = Ki, where we get the Barycentric refinement operator A, the eigenvec¬ 
tor of to the eigenvalue 1 gives an invariant for Barycentric refinement: 
it is the Euler characteristic vector (1,—1,1,—1,1,—1,...). When restricted 
to d-graphs, the vector ( 0 , 0 ,..., 1 ) with 1 at the d + I’th entry is an eigen¬ 
vector to the eigenvalue (d-l- 1)!. It leads to a counting invariant in the limit 
which is called volume. The other eigenvectors lead to similar limits which 
must be linear combinations of valuations in Riemannian geometry |56ll32j . 
We see that the even Barycentric invariants X 2 k{G) = f 2 k{A) ■ v{G) 
are constant zero on d-graphs, where fk{A) are the eigenvectors of A'^. It 
will imply for example that for any compact smooth 4-manifold contain¬ 
ing V vertices, e edges, / triangles, g tetrahedra, and h pentatopes satisfies 
22e -|- 40^ = 33/ -|- 45/i. In the continuum, the more general operator Ah 
leads to a linear map on valuations. These maps will help to investigate the 
connection between the discrete and continuum. 

7) In the same way as for Riemannian manifolds, where one studies the 
eigenfunctions fk of the Laplacians, the eigenfunctions of the Kirchhoff 
Laplacians will play an important role for understanding the limiting den¬ 
sity of states. The Chladni figures are the level sets fk = 0. They can be 
defined for d-graphs. We can look at / = 0. It is by a discrete Sard result 
[H] a (d — l)-graph as long as 0 is not a value taken by /. If 0 is a value 
taken we just can plot / = e for 0 < e small. The point is that we can so 
define nice (d — 1) graphs associated to fk- The topology of these Chladni 
graphs depends very much on the energy as in the continuum and as in the 
continuum, are not well understood yet. 

8 ) The even eigenvectors V 2 k of A^ appear to have the property that X 2 k{G) = 

V 2 k'v{G) = 0 for any 2d-graph. If true, this leads to integral geometric invari¬ 
ants in the limit. For any graph, define X 2 k{G) = limm->.oo V 2 kv{Gm)/■ 
The limit is trivial as the right hand side is constant. For d-graphs, we see 
X 2 k{G) = 0. We have tried random versions of 5^, 5^ x 5^ x T^, x 

, 5®. The limiting integral theoretic invariants for d-manifolds would be 
zero for differentiable manifolds. Assume we would find a graph G with 
uniform dimension d which is not a d-graph but which is homeomorphic 
to a d-graph H. Then, since Barycentric refinements preserve the homeo- 
morphism relation, the Barycentric limit M of G is a topological manifold 
homeomorphic to the Barycentric limit N of H. While M, N are topolog¬ 
ical manifolds which are homeomorphic, they can not be diffeomorphic, as 
integral geometric integer valued invariants are diffeomorphism invariants. 
A basis for valuations can be obtained if M is embedded in a projective 
n-sphere. We can compute the expectation of the /c-volume of a random 
m-planes with M using a natural probability measure obtained from Haar 
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measure on SO{n) acting on /c-planes. This leads to d + 1 invariants, for 
which some linear combinations is the Euler characteristic. The graph the¬ 
oretical invariants obtained from A could produce invariants allowing to 
distinguish homeomorphic but not diffeomorphic manifolds. 


5. About the literature 

Barycentric rehnement are of central importance in topology. Dieudonne 
|16j writes of ’’the three essential innovations that launched combinatorial 
topology: simplicial subdivisions by the barycentric method, the use of dual 
triangulation and, finally, the use of incidence matrices and of their reduc¬ 
tion. ” In algebraic topology, it is used for proving the excision theorem or 
the simplicial approximation theorem |25L I55| . In topology, Barycen¬ 
tric subdivision is primarily used in an Euclidean setting for subdividing 
complex polytops or CW complexes. For abstract simplicial complexes, it is 
related to flag complexes even so there are various inequivalent definitions 
of what simplicial subdivision or simplex graph or Barycentric subdivisions 
are. In graph theory, subdivisions are considered for simplicial complexes 
which are special hypergraphs in topological graph theory [23]. Indeed, 
most graph theory treats graphs as one-dimensional simplicial complexes, 
ignoring the Whitney complex of all complete subgraphs. Subdivisions 
classically considered for graphs only agree with the definition used here if 
the graph has no triangles. Two graphs are classically homeomorphic, if 
they have isomorphic subdivisions ignoring triangles (see e.g. |ll[30l[2]). In 
the context of maps which are finite cell complexes whose topological space 
is a surface S, Barycentric subdivisions are considered for this cell complex 
|22| but not for the graph. In the context of convex polytops, Barycentric 
subdivision appear for the order complex of a polytop which is an abstract 
simplicial complex |48| . 

Resources on the spectral theory of graph are in PESKUdalEIlEII. it 
parallels to a great deal the corresponding theory for Riemannian manifolds 

[TlEllEI. 

The definition of d-spheres and homotopy are both due to Evako. Having 
developed the sphere notion independently in miiiai we realized in |l2| the 
earlier definition of Ivashchenko=Evako |29l EH] EZl EB] . The definition of 
these Evako spheres is based on Ivashchenko homotopy which is homotopy 
notion inspired by Whitehead [^3] but defined for graphs. 

The Barycentric refinements for d = 2 are studied in |5], where the limit of 
Barycentric refinement has a dual called the hexacarpet. Figure 7) in that 
paper shows the eigenvalue counting function in which gaps in the eigen¬ 
values are shown similarly as in the spectral function Fq{x). The work [5] 
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must therefore be credited for the experimental discovery of the gap. 

Papers like [13 ES] deal with the geometry of tessellations of a triangle, 
which defines a random walk leading to a dense subgroup of SL{2, R) which 
defines a Lyapunov exponent. While the focus of those papers is different, 
there might be relations. 

For illustrations, more motivation and background, see also our first write¬ 
up go]. 


6. Figures and Code 

Here are the Mathematica routines which produced the graphs: 


TopDim = 2; 

Cl [ s _ , k_ ]: = Module [ { n , t ,m, u , q , V—V ertexList [s] ,W^EdgeL ist [ s ] , 1 } , 
n—Length[V]; rr^Length [W] ; u—Subsets [V, { k , k }] ; q—Length[u]; 1^{}; 
VWTable[{W[[j . 1 ] ] .W[ [ j . 2 ] ] } . { j .m} ] ; I f [ k==l. l=Table [ { V [ [ j ] ] } .{j ,n}] , 
If [ k==2, l^V, Do [ t=Subgraph [s,u[[j]]]; If[ Length [EdgeList [t]]= = 
Binomial [ k , 2 ] , 1—Append [1 ,VertexList[t]]] , {j ,q}]]]; 1]; 

Ring[s_ ,a_]: = Module [{v,n,m,u,X} , v=Ver texL ist [ s ] ; n=Length [ v ] ; 
u—Table [Cl[s,k],{k, TopDim +1}] / . Table [ k—>a [[k]] ,{k,n}] ;rr^Length [ u ] ; 

X^um [Sum[ Product [u[[k,l,m]] , 

{m, Length [u [ [k , 1 ] ] ] }] , { 1 .Length [u [ [k]]]}] ,{k,m}]]; 

GR[ f - ]: = Module [ { s = { } } ,Do[Do [ If [ Denominator [f[[k]]/f[[l]]] = =l k! = 1 , 
s—Append [ s , k—>1 ] ] , { k , Length [ f ] } ] , { 1 , Length [ f ] } ] ; 

UndirectedGraph [ Graph [ s ] ] ] ; 

GraphProduct [ s 1 _ , s2 _ ]: — Module [{f,g,i ,fc ,tc}, 
f c=FromCharacterCode ; t c=ToCharacterCode ; 

i [ 1 _ , n_]: — Table [ fc [Join [ tc [ 1 ] , IntegerDigits [k] + 48]] ,{k,n}] ; 
f—Ring [si , i [ ” a” , Length [VertexList [si ] ] ] ] ; 

g—Ring [ s2 , i [ ” b” , Length [VertexList [s2 ] ] ] ] ; GR[ Expand [ f *g ] ] ] ; 

NewGraph [ s_ ] : = GraphProduct [ s , CompleteGraph [ 1 ] ] ; 

Bary[s_ ,n_]: = Last[NestList[ NewGraph , s , n ] ] ; 

Laplace [s_ ,n_]: = Normal [KirchhoffMatrix[Bary[s ,n]]] ; 

F[s_ ,n_]: = Module [ { u ,m} , u=Sort [Eigenvalues[1.0*Laplace [s ,n]]] ; 

rr^Length [u ] ; Plot [u [ [ Floor [x*m] ]] ,{x,0 ,1}]]; 

TopDim— 1; Show [ Table [F[CycleGraph[4] ,k] ,{k,3}] , Plot Range — >{0,4}] 
TopDim —2; Show [ Table [F[CompleteGraph[3] ,k] ,{k,3}] , Plot Range —>{0,10}] 


And here is the recursive computation of the matrix A which gives the 
clique data Av of the Barycentric refinement of Gi if v is the clique vector 
of G. The computation produces the clique data of the boundary of {Kn)i 
which allows to compute the number of interior /c-simplices producing the 
off diagonal matrix elements not in the top row. 


BarycentricOperator [m_]: = Module[{ } , 

b [ A_]: — Module [ { n—Length [A] , c } , c—A. Table [ Binomial [n+l,k] ,{k,n}]; 

Delete [ Prepend [c,l] ,n+l]]; 

T [ A_] := Append [ Transpose [Append [ Transpose [A] ,b[A]]] , 

Append [ Table [0 , { Length [A] } ] , ( Length [A] + 1 ) ! ] ] ; 

Last [NestList [T,{{1}} ,m] ] ] ; BarycentricOperator [7] 
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Figure 1. The functions for G = converge to 
the limiting function 4sin^(7rx/2). The functions for 

G = converge to a limiting function which appears to 
have jumps corresponding to gaps in the spectrum. We have 
here only established that the universal limit F exists, but 
already for d = 2 do not know about the nature of the limit. 
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The procedure produces finite dimensional versions of this matrix which 


matter when 


pie, we have 


ooking for interesting quantities on Qd- For d = 2 for exam- 
111 ^ 

whose transpose has the eigenvectors [1, —1,1]^ 


(Euler characteristic), the eigenvectors [0, —2, 3]^, a well known invariant for 


2-graphs which is proportional to the boundary for 2-grap 


IS with boundary 


as well as [0, 0,1]^ which is area. For d = 3, where A = 

the eigenvectors besides the Euler characteristic vector [1, —1,1, —1]^ to the 
eigenvalue 1 and volume [0, 0, 0,1]^ to the eigenvalue 24, there is [0, 22, —33,40]^ 
with eigenvalue 2 and [0, 0, 0 — 1,2]^ with eigenvalue 6. The later gives an 
invariant which is zero for 3-graphs as in general [0,... 0, —2, d -|- 1] is an 
invariant for d-graphs. We see experimentally that for any d-graph,the 
eigenfunction of to each eigenvalue (2A;)! is perpendicular to 
the clique vector v. For example, there is a discrete x 5^ with 
clique vector v = [1908, 26520,87020,104010,41604] which is perpendicular 
to V 4 = [0, —22, 33, —40,45]. The invariant does not change under edge re¬ 
finement modifications which are homotopies preserving d-graphs. The 
invariants X 2 k{G) = v{G)-V 2 k remain zero under such homotopies. They are 
currently under investigation. It looks promising that integral theoretical 
methods like |37l [39] for generalized curvatures allow to prove the invariants 
to be zero. It would be useful also to know whether any two homeomorphic 
d-graphs have a common refinement when using Barycentric or edge refine¬ 
ments. This looks more accessible than related questions for triangulations 
as d-graphs can be dealt with recursively and refinements of unit spheres 
carry to refinements of the entire graph. 
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